Field induced linear magneto-elastic coupling in multiferroic TbMnOa 
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We have used in-field neutron and X-ray single crystal diffraction to measure the incommensu- 
rability S of the crystal and magnetic structure of multiferroic TbMnOa. We show that the flop 
in the electric polarization at the critical field He, for field H along the a— and b— axis coincides 
with a fst order transition to a commensurate phase with propagation vector k = (0, i, 0). In-field 
X-ray diffraction measurements show that the quadratic magneto-elastic coupling breaks down with 
applied field as shown by the observation of the 1st harmonic lattice reflections above and below 
He- This indicates that magnetic field induces a linear magneto-elastic coupling. We argue that 
the commensurate phase can be described by an ordering of Mn-O-Mn bond angles. 
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Control of the spontaneous ferroelectric polarization 
(P s ) with an external magnetic field (if) in a material 
opens the opportunity for new types of magneto-electric 
devices. The realization of such devices is based on multi- 
ferroic materials in which magnetism and ferroelectricity 
are strongly coupled. While available multiferroics are 
limited, it has been shown that frustrated spin materials 
may offer a new class of enhanced multiferroics 0, ■ 
In one of these new materials, TbMnC>3, we find that mul- 
tiferroic behavior arises as a consequence of the release of 
frustration with if. Here ferroelectricity arises below the 
Neel temperature (Tjy), from a coupling to the lattice 
of an incommensurate (1C) modulation of the magnetic 
structure (Fig. ^i) that is caused from frustration in 
the ordering Mn d orbitals [ij, |3j ■ In this letter we show 
that magnetic field releases this frustration, inducing a 
linear magneto-elastic coupling, so that ferroelectricity is 
no longer a secondary effect. The linear magneto-elastic 
coupling drives a magneto-structural transition from an 
IC phase which has P s along the c— axis to a commensu- 
rate (C) phase with P s along the a-axis. 

In TbMn0 3 (orthorhombic, a=5.302A, 6=5.8571, 
c=7.402A) when a magnetic field is applied along the 
6-axis (if 1 1 6) at 2K, parallel to the direction of the IC 
magnetic modulation (see Fig. |TJi), the electric polar- 
ization of the lattice flops from along the c— axis (P||c) 
to along the a— axis (P\\a) at the critical field if^-4.5T 
(Fig-Eld)). When field is applied along the a-axis (if \\a), 
perpendicular to the magnetic modulation, a similar flop 
is found but at a higher critical field, ffg~9T (Fig.^Ja)). 
Recently there have been a number of examples of novel 
magneto-elastic coupling in complex multiferroic oxides 
such as Tb2Mn205 which exhibits a reversable polariza- 
tion switch with applied field 0| , and hexagonal HoMn03 
where one magnetic phase is selected over another by ap- 
plying an electric field y|. However, TbMn03 is unique 
as it is the only known example of a material that ex- 
hibits a field induced flop of the polarization. 

In TbMn03 the staggered ordering of Mn 3+ 



3d,3 X 2_ r 2 /3d 3y 2_ r 2 orbitals as found in LaMnC>3 is frus- 
trated partly due to the small ionic size of Tb 3+ 0|. 
This leads to an IC spin ordering which drives a fer- 
roelectric lattice modulation 0, Q . The wavevector for 
the modulation of the Mn magnetic moments is K^ n = 
(0, k±5™ n , 0), with incommensurability (5* fn ~0.29 at T^v 
(1st harmonic) Accompanying the magnetic or- 
dering there is a lattice modulation with <5/=2<5 m (2nd 
harmonic) 0, 0, consistent with a quadratic magneto- 
elastic coupling between the lattice and the spin density 
wave 0,^3- Below T-jy b =7K, sinusoidal Tb spin order 
is also found (k™ = (0, k±5™, 0), 5™=0.42) 8] and is 
expected to produce a similar coupling to the lattice with 
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To uncover the nature of the coupling between P and 
if we have used in-field neutron and synchrotron X-ray 
single crystal diffraction to measure the field response 
of S with field applied parallel to the a— and 6— axis 
(if || a and if||6), up to 14T. Single crystals were ob- 
tained by re-crystalizing a ceramic rod of TbMnC>3 un- 
der Ar atmosphere using an optical floating zone furnace. 
Small pieces of our crystal were characterized using mag- 
netization and specific heat with good agreement with 
published measurements 0. Synchrotron X-ray diffrac- 
tion measurements were carried out on beamline X21 
(Ei =9.5 keV) at the NSLS, Brookhaven National Labo- 
ratory, using a vertical 13T superconducting magnet with 
the sample alighed with if||a, and at the BW5 beamline 
at HASYLAB (JSj =100.5 keV), using a horizontal 10T 
superconducting magnet and the sample aligned with 
if || b. Neutron diffraction experiments were made on a 
single crystal of TbMn03, (19mm in length and 8mm 
in diameter) at the BENSC facility of the Hahn-Meitner 
Institut, Berlin. For measurements with if||a we used 
the FLEX cold triple-axis spectrometer with collimation 
of 40'-40'-40', fci=1.5A -1 , and a cooled Be filter on the 
scattered beam, while a magnetic field was applied using 
the vertical 14. 5T superconducting magnet VM1. For 
measurements with if || b we used the 2-axis diffractomc- 
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ter E4 (A = 2.2 A) with 60'-40' collimation while field 
was applied using the 4T horizontal field magnet HM1. 
In all cases diffraction measurements were made in the 
frc-plane. Our neutron diffraction measurements are sen- 
sitive to magnetic ordering whereas X-ray diffraction ex- 
periments probe the lattice modulation. 

We applied field parallel to the a— axis and we mea- 
sured characteristic reflections using neutron diffraction 
by cooling our crystal from 50 to 2K in fields of 0,2,4,6,10 
and 14T. The variation of S for Mn and Tb spin ordering 
as a function of field at 2K, compiled from these mea- 
surements, is shown in Fig. Elk), while the variation of 
the intensity of the same reflections is shown in Fig. EI C )- 

In Fig. Etb) we show neutron diffraction scans through 
the 1st harmonic reflection (0,i5^f ra ,l) as a function of 
H\\a. The intensity of this reflection is enhanced with 
magnetic field from 2 to 4T while the magnitude of <5*f™ 
remains unchanged for fields up to 8T (FigEtb,c)). How- 
ever, for H >H C ~9T we find that S of the 1st harmonic 
reflection has changed to a C value of 5=\ and is ac- 
companied by a 2nd harmonic reflection with 28=\ (Fig. 
Efb)). This transition at H c coincides with the polariza- 
tion flop from P\\c to P\\a as shown in Fig. Et a )- The 
implication of this observation is that the polarization 
flop arises because the spontaneous polarization of the 
C-phase is along the a— axis, different to that of the IC- 
phase. 

In Fig. Eta) we show X-ray diffraction scans as a func- 
tion of H\\a at T=10K for the 1st harmonic (0,3-5 Mn ,8) 
reflection. If the magneto-elastic coupling is quadratic, 
then only magnetic scattering is allowed for this reflec- 
tion, and indeed we do not observe a signal at this posi- 
tion for H\\a<2T. However, for H\\a>2T a peak becomes 
visible and above Hq it locks into the C value of 5=4. 
The observation of a lattice reflection at the wavevec- 
tor of the 1st harmonic indicates that for H\\a a linear 
magneto-elastic coupling is induced above 2T which is 
maintained above He into the C-phase. The variation of 
the intensity of this reflection, shown in Fig. Et c ) is com- 
plex, but its changes are associated with the metamag- 
netic transition of Tb-spins (see below) and the increase 
in the intensity of the 1st harmonic reflection (0,6 m , 1) 
between 2 and 4T (see Fig. Etc)). 

The behavior of the magnetic reflections due to Tb or- 
dering is significantly different from that for Mn. The in- 
tensity of the 1st harmonic reflection (0,5™, 1) decreases 
significantly between and 2T (Fig. Etc)), indicative of 
a metamagnetic transition reported previously fl. This 
transition in the Tb magnetic sub-lattice is coincident 
with an anomaly in P\\c as well as with a sharp increase 
in the magnetization reported elsewhere Q, How- 
ever, above H c we find no evidence of sinusoidal Tb spin 
ordering, suggesting that Tb spins have also ordered com- 
mensurately. Rather we find that with H\\a Tb spins are 
polarized ferromagneticaly as indicated by an increase in 
the intensity of Bragg reflections. 



The polarization flop when the field is applied along the 
b— axis in the direction of the magnetic modulation, oc- 
curs at a lower critical field of H C ^4T (Fig. Eld)). Here 
we applied a magnetic field parallel to the 6— axis after 
the sample was cooled to 2K in zero field and characteris- 
tic superlattice reflections, such as the (0, 2S Mn , 3), were 
measured using synchrotron X-ray diffraction (see Fig. 
4). When a magnetic field is applied, the incommensu- 
rability of the 2nd harmonic reflection remains invariant 
(Fig. Et e )) but its intensity rapidly decreases with field 
from to 2T and remains low up to ~4T ( Fig. EfO)- I n ~ 
creasing field clearly destabilizes the lattice modulation 
and leads to a discontinuous transition to the C-phase 
with 28=1 at #c= 4 - 5T ( Fi §- Et a ) andEte,f)). We find 
C and IC reflections coexisting over the field range of 4 to 
5.5T, while for H >5.5T only the 26=^ superlattice re- 
flection is visible (Fig. Eta)). Coincident to this magneto- 
structural transition is the polarization flop from P\\c to 
P\\a at about 4.5T as shown in Fig.Etd)- 

As with the previous field configuration, here also 
our X-ray measurements revealed superlattice reflections 
with S=\ above H b c (Fig. \gb) and Et e , f ))- Th c lin- 
ear increase in the intensity of this reflection with H\\b 
above H h c is particularly striking, as it suggests a linear 
magneto-electric coupling in the C-phase. The variation 
of intensity with field of the 2nd harmonic reflection is 
different in that it increases rapidly above H c and shows 
a constant variation with H above ~6T (Fig. Elf))- 

For this field configuration we find using neutron 
diffraction that at 2K Tb spins undergo an additional 
discontinuous transition to a phase with S™—^ at 
-ff ||&=1.25T, as shown in Fig.Et e )- From symmetry anal- 
ysis Tb-spins for this phase are both sign and amplitude 
modulated. This transition is reflected in an anomaly in 
P\\c at the same field value (Fig.Etd)), and a significant 
increase in the magnetization at this temperature is also 
reported^. Above H c we find no evidence of incommen- 
surate ordering of Tb-spins, but the significant increase 
of the magnetization at this field indicates that Tb spins 
order commensurately there 0. 

In zero field TbMn03 is an improper ferroelectric 
where the lattice polarization arises from a quadratic 
magneto-elastic coupling to a spin density wave as indi- 
cated by the observation of the 2nd harmonic structural 
reflections 0, . Here we demonstrate that this picture 
of improper ferroelectricity in TbMn03 breaks down with 
H even below He for the case of H\\a and most likely 
for H\\b. The identical periodicity of the magnetic and 
crystal structure indicates that H switches on a linear 
magneto-elastic coupling which is maintained into the C- 
phase. This would suggest that the IC-C phase transi- 
tion we report on and consequently the polarization flop 
at He are driven by field induced linear magneto-elastic 
coupling. Further our measurements show that although 
the application of field ultimately results in a C-phase, 
the effect of H on Tb-spins is dependent on the field 
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direction. For H\\b the transition to the C-phase is cou- 
pled to a metamagnetic transition of Tb-spins Q , while 
for H\\a the Tb-spins order commensurately at field be- 
low Hq. How exactly Tb-spins couple to the transition 
of the high field C-phase is yet to be resolved. 

The magneto- electric coupling is expected to change 
in a similar way as the magneto-elastic coupling. Indeed 
the behavior of TbMnC>3 with H is reminiscent to that of 
BiFeOa which exhibits a switch to linear magneto-electric 
coupling at ~20T^^- The coupling between P and H is 
expressed by Pi — P S i + aijHj + ^PijkHjH^, where P s 
is the spontaneous polarization, cnj is the tensor of the 
magneto-electric susceptibility, ftijk is a tensor describing 
the coupling of an incommensurate magnetic structure 
to the lattice 01 • I* 1 an analogy to BiFeOa, we suggest 
that as field destroys the incommensurate spin ordering 
(Pijk — *0), a linear magneto-electric effect would be ex- 
pected to arise Indeed a linear behavior in P vs. 
H is evident for H\\b (Fig. 2(d)). 

The commensurate propagation vector of k — (0, i, 0) 
uniquely gives a spin structure that is only sign modu- 
lated as depicted in Fig. QJb) . It is of interest to consider 
the effect of this magnetic ordering on the frustrated or- 
bital ordering. The crystal structure of TbMn0 3 (H=0T, 
T >Tjv) shows an ordering of long and short bonds as 
found in LaMn03, but with a low angle (ft of 144° be- 
tween Mn06 octahedra (Mn-O-Mn), compared to 165° 
found in LaMnOs [l4| . A simple application of Goode- 
nough's rules would suggest a layered antiferromagnetic 
(AF) ordering, inconsistent with the present data[l5| . 
We propose here that the structural distortion that is 
coupled to the till spin structure arises from the or- 
dering of (ft. Small values of (ft would result in AF in- 
teractions between adjacent Mn-spins as the influence of 
su per- exchange via fog-orbitals would be more significant 
[3lll5l| . On the other hand a relative larger value of (ft will 
be less influenced from fog interactions and would prop- 
agate a ferromagnetic coupling between Mn-spins as it 
does in LaMnC>3. Indeed the layered AF structure is sta- 
ble to values of (ft as low as 148°, found in SmMnO-3 0. 
Using the magnetic ordering of the high field C-phase as 
a constraint and displacing O-atoms in the para-electric 
structure so that (ft should increase (decrease) for ferro- 
magnetic (antiferromagnetic) coupling between adjacent 
Mn-spins, one arrives at a structural model show in Fig. 
^d. The result is a structure were ferromagnetically cou- 
pled MnOe octahedra undergo simple rotations (ferro- 
electrically inactive), while octahedra that are frustrated 
by antiferromagnetic coupling in the a6-plane show a 
scissor-like distortion (ferroelectrically active). The com- 
bined rotations and scissor-like distortions in the struc- 
ture leads to anti-ferroelectric displacements for the O- 
atoms along the b— axis but an overall ferroelectric dis- 
placement along the a— axis in agreement with polariza- 
tion measurements. While this simple model predicts the 
correct polarization for the C-phase on the basis of the 



magnetic structure, the periodicity of the structural dis- 
tortion suggests displacements of the Mn ions that can 
not be accurately predicted in the absence of high field 
crystallographic data. 

In this letter we have demonstrated that the polariza- 
tion flop in TbMn03 arises from a transition to a C-phase 
with propagation vector n — (0, 4,0). The quadratic 
magneto-elastic coupling breaks down with field and 
the high field phase exhibits linear magneto-elastic and 
magneto-electric behavior. We argue that field has the 
effect of releasing the apparent frustration in the orbital 
ordering in TbMnC>3 and leads to a C-phase by an or- 
dering of Mn-O-Mn bond angles. 
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FIG. 1: Crystal and magnetic structure of TbMnOs in the 
ab— plane. Long (darker) and short (lighter) Mn-O bond or- 
dering arises from staggered 3d 3;[ ,2_ r .2 /3d 3y 2_ r 2 orbital order- 
ing found above Tat. (a) A Model of the IC magnetic struc- 
ture of TbMnC>3 below Tjv Q. Arrows through Mn atoms 
indicate the direction of Mn-spins. (b) A model of the Mn- 
spin structure for TbMnOs for H>i/c- Here arrows from 
O-atoms indicate the direction of displacement from the av- 
erage structure as to increase (decrease) the angle <j) on the 
basis of ferromagnetic (antiferromagnetic) coupling between 
adjacent Mn-ions. Large arrows on ferroelectrically active oc- 
tahedra indicate the direction of the predicted polarization 
from these displacements. 



FIG. 3: (a) X-ray diffraction scans of the 1st harmonic (0,3- 
K,8) reflection as a function of field with H\\a at T=10K, 
measured on X21. (b) Neutron diffraction scans through the 
1st harmonic (0,<5^f™,l) reflection as a function of H\\a. 



FIG. 2: (a,d) Electric polarization of TbMnOs parallel to the 
a- (P\\a) and c-axis (P||c) with H\\a and H\\b, respectively. 
These data were compiled from field cooled measurements 
reported by Kimura et aZ. (b) Variation of S as a function 
of H\\a, for the 1st harmonic reflections (0,6^ n , 1), (0,5™ ,1) 
measured using neutron diffraction (filled blue and red circles 
respectively) and (0, 3—5, 8) measured with X-ray diffraction 
(filled green squares). For H >Hq the C reflections (0, j, 1) 
and (0, ~ j 1) (open blue circles) were measured using neutron 
diffraction while the (0, 2|,8) reflection with X-rays (open 
green squares), (c) Variation of the integrated intensity as a 
function H\\a, for the same reflections as in panel (b). For 
the reflection due to Tb ordering the logarithmic scale on the 
right axis is used, (e) Variation of 5 as a function H\\b, for 
the 1st harmonic reflections (0,<J^ ,1) and C (0, |, 1) due to 
Tb spin ordering (filled and unfilled red circles respectively). 
On the same panel we show the results of X-ray diffraction 
measurements of the 2nd harmonic Mn reflection (0, 25 , 3) 
(filled blue circles) and the 1st and 2nd harmonic reflections 
(0, |, 3) and (0, |, 3) (open blue circles), (f). Variation of the 
integrated intensity as a function of H\\b, for the reflections 
shown in panel (e). For the reflection due to Tb ordering 
the logarithmic scale on the right axis is used. Integrated 
intensity in panels (c) and (d) is shown in arbitrary units. 



FIG. 4: (a) Scans of the 2nd harmonic IC lattice reflection at 
(0, 2(5 M ",3) with increasing field at 2.6K measured on BW5. 
The C reflection at (0,i,3) is shown for H >Hq. (b) Similarly, 
scans of the 1st harmonic (0,|,4) reflection of the C-phase 
with increasing field at 2.6K for H >Hq are shown. 
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